NASA TN D-416 


N E2 70990 

NASA TN D-416 



TECHNICAL NOTE 

D-416 


ANALYTICAL INVESTIGATION OF THE DYNAMIC BEHAVIOR 

OF A NONLIFTING MANNED REENTRY VEHICLE 

By Jacob H. Lichtenstein 

Langley Research Center 
Langley Field, Va. 


i 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON September 1960 




NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


TECHNICAL NOTE D-Ll6 

ANALYTICAL INVESTIGATION OF IHE DYNAMIC BEHAVIOR 
OF A NONLIFTING MANNED REENTRY VEHICLE 
By Jacob H. Lichtenstein 


SUMMARY 


An analytic investigation was made of the dynamic behavior of a 
nonlifting manned reentry vehicle as it descended through the atmosphere . 
The investigation included the effects of variations in the aerodynamic 
stability derivatives, the spin rate, reentry angle, and velocity. The 
effect of geostrophic winds and of employing a drogue parachute for 
stability purposes were also investigated. 

It was found that for the portion of the flight above a Mach number 
of 1 a moderate amount of negative damping could be tolerated but below 
a Mach number of 1 good damping is necessary. The low-speed stability 
could be improved by employing a drogue parachute. The effectiveness 
of the drogue parachute was increased when attached around the periphery 
of the rear of the vehicle rather than at the center. Neither moderate 
amounts of spin or the geostrophic winds had appreciable effects on the 
stability of the vehicle. The geostrophic winds and the reentry angle 
or velocity all showed important effects on the range covered by the 
reentry flight path. 


INTRODUCTION 


There is, at present, considerable interest in the dy nam ics of 
bodies entering the earth's atmosphere. This interest stems from the 
concern about the behavior of both satellite (manned or unmanned) and 
ballistic missiles as they descend through the atmosphere. As a result 
of this interest numerous papers have been published by the NASA and 
others, some of which discuss the trajectory of reentry bodies (for 
instance, refs. 1 to A) and some of which discuss the dynamic stability 
of the reentry body as well (for instance, refs. 5 to 12). 

The present paper reports the results of an analytic investigation 
into the dynamic behavior of a representative nonlifting manned satellite 
as it descends through the atmosphere. This investigation encompassed 
the effects of changes of the aerodynamic derivatives and rate of spin 
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on the stability of the vehicle. Computations were made to assess the 
advantages of using a drogue parachute to enhance the subsonic stability 
of an otherwise marginal or unstable body. In this connection the rela- 
tive merits of attaching the parachute towline at one central point or 
at several points around the periphery at the rear of the body were also 
investigated. In addition, the effects of the earth's spin and of varia- 
tions in the angle of entry and initial velocity were computed. The 
computations for this program were made on an IBM 70^ electronic data 
processing machine. 


SYMBOLS 


The axes system used in the present program is shown in figure 1. 
A maximum cross-sectional area, sq ft 

a radial distance from center line of body to parachute 

attachment point on body, ft 

a used as step input at some specified value of 
with the sign of a varying as the sign of 

(see appendix A) 

a used as step input at some specified value of 
with the sign of a varying as the sign of 

(see appendix A) 
drag coefficient 


value of 

v - V v 


Vr 

v - V„ 


'R 




value of 

v - v z 

V R 

W - V-7 


'R 



Pitching moment 
QAd 


variation of pitching- moment coefficient with angle of 

5c 


attack in XZ-plane, 



per radian 




damping-in-pitch coefficient, 


-q o\ oo f 



-q on cd t-f 
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p _ Yawing moment 
n " QAd 


variation of yawing-moment coefficient with angle of attack in 

dC n 

XY -plane , — -, per radian 

C^Oy 


damping-in-yaw coefficient, 


Axial force 


Axial force of drogue parachute 


Side force 


variation of side-force coefficient with nondimens ional 

c>C Y 

yawing velocity, 

a i: 


C Y rate of change of side-force coefficient with angle of 

a SC Y 

attack in XY-plane, — — , per radian 

oay 

P Normal force 
^ QA 

C z rate of change of normal-force coefficient with angle of 

a 

attack in XZ-plane, - — , per radian 

“Z 


variation of normal-force coefficient with nondimens ional 


pitching velocity, 
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a 

F 

F 

r g 

F X> F Y 

g 

h 


I Y^ I Z 


K 

l 


M 

m 


maximum body diameter, ft 
force, lb 

force due to gravity, lb 

forces along X- and Y-axes, respectively, lb 

acceleration due to gravity at surface of earth, ft/sec^ 
altitude above surface of earth, R - R e , ft 

p 

moment of inertia, slug -ft 

p 

moment of inertia about X-axis, slug-ft 

moment of inertia about Y- and Z-axes (same value), 
respectively, slug-ft^ 

,2 

damping factor, -2C X + 

distance rearward from body center of gravity to drogue- 
parachute attachment point on body, ft (see sketch in 
appendix B) 

moment about center of gravity 

moments about X-, Y-, and Z-axes, respectively, ft-lb 
mass, slugs 


%a 


n 




Q 


Mach number, - 

Speed of sound 

acceleration along flight path, g units 

angular velocities about X-, Y-, and Z-axes, respectively, 
radians /sec 

dynamic pressure, ipVp , lb/sq ft 


PoO^O ^ 



^ oo vo r~- 


5 


qd_ 

2V R 

R 


Re 

rd 

2Vr 

u,v,w 


V 

V R 

V X> V Y* V Z 

X,Y,Z 


Xi,Yi,Zi 

x,y,z 


CXy 

“z 


nondimens ional pitching-velocity parameter referred to 
diameter 

radial distance from origin of the Xf-, Yf-, and Z-^-axes 
system to origin of X-, Y-, and Z-axes system, ft 

radius of earth, taken as 21 x 10^, ft 

nondimens ional yawing -velocity parameter referred to 
diameter 

components of velocity along X-, Y-, and Z-axes, respectively, 
ft/sec 

resultant inertial velocity of vehicle, tfu^'+'v 2 + w 2, ft /sec 

resultant velocity of vehicle with respect to air, ft/sec 

components of geostrophic wind along X-, Y-, and Z-axes, 
respectively, ft /sec 

reference axes with origin at center of mass of body with 
X-axis alined along line of symmetry 

inertial reference axes with origin at center of earth 

distance from origin along X^-, Y-^-, and Z^-axes to center 
of gravity of vehicle 

total angle of attack, deg 

angle of attack in XY-plane, deg 

angle of attack in XZ-plane, deg 


7 flight-path angle, deg 

Euler angles defined in figure 1, radians 


0 angle between X^Z— plane and location vector to vehicle 

(latitude ) 

density of air (standard atmosphere tables for density 
variation with altitude, from Rocket Panel Proposal of 
April 1955); slugs /cu ft 


P 
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o 


radius of gyration 



sq ft 


\|r angle in XpZi-plane between XjYj -plane and plane containing 

vehicle (west longitude) 

ft rate of spin of the earth, taken as 0.000073 radian/sec 

Subscripts : 

max maximum 

min minimum 

o indicates initial value 

p indicates drogue parachute 

A dot over a symbol indicates differentiation with respect to time. 


CALCULATIONS 


Computations were made of the motion of a body representative of a 
manned satellite upon its reentry into and descent through the earth's 
atmosphere. The computations were made on an IBM 704 electronic data 
processing machine and the conditions for the computations are given in 
table I. The dimensional and mass characteristics of the vehicle are 
given in table II. The general reentry configuration considered for the 
study is shown in figure 2 and the methods of drogue -parachute attach- 
ment are shown in figure 3 • For the two-point attachment with 
a/d = 0.3125 the parachute lines are attached at the periphery of the 
rear face. For the two-point attachment with a/d = 0.625 the parachute 
lines are attached to arms extending outwarc from the surface at the 
rear. The equations of motion used in these calculations are given in 
appendix A and the development of the terms used to simulate the drogue 
parachute are given in appendix B. 


For a computation of this type where the frequency of oscillation 
varies through very wide limits (Frequency ^/q”), a constant computing 
interval which would insure repeatability at the high frequencies would 
be very wasteful at the lower frequencies. Therefore, a program in 
which the number of points per oscillation cycle was prescribed was used 
rather than a fixed time interval. It was found that computing between 
20 and kO points per cycle generally was sufficient for adequate 
repeatability. 


1 oo vo 
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The following general pattern was used in making the computations: 
From reentry at an altitude of 380,000 feet down to an altitude where 
%[a = the aerodynamic derivatives were considered constant and the 
drogue -parachute terms were considered zero} below a Mach number of 1 
some aerodynamic derivatives were changed to values that were more 
representative of low-speed values, and the drogue -parachute terms were 
included. Thus, it was possible to compute for the low-speed part of 
the run the various effects of aerodynamic derivatives and drogue- 
parachute configuration without the necessity of recomputing any of the 
high-speed part of the run. The mass characteristics and stability 
derivatives used were the most representative of the desired configura- 
tions that were available at the time the computations were made. Some 
computations were made with the mass characteristics shown for condi- 
tion 1 (see table II) which were early estimates. These values were 
changed to those listed for condition 2 when these more representative 
data became available; however, it was not deemed worthwhile to recompute 
the earlier runs . 


RESULTS AND DISCUSSION 


Presentation of Results 

The results of this investigation are presented in figures 4 to 14 
according to the arrangement shown in table I. The data divide con- 
veniently into two parts: the high-speed portion, from reentry down to 

the altitude where Nj^ a = 1, and the low-speed portion, below N^ a = 1 

where the drogue parachute generally was employed. The data for the 
high-speed portion are presented in figures 4 to 9. The results of the 
low-speed portion of the trajectory, including the effects of employing 
a drogue parachute for stability, are shown in figures 10 to 14. Gen- 
erally, these computations were made starting with the pbint at Nj^ a = 1 

for the basic condition (run 1, table I). 


Basic Condition 

The data for the basic condition are given in figure 4. This con- 
dition is listed as run 1 in table I and the aerodynamic derivatives that 
were used are C m = 0.37, C™ = -0.20, and C z = -0.28 for Mach num- 

bers greater than 1 and Cn^ = -1.00, = - 0 . 05 , and C z ^ = 0.60 for 

Mach numbers less than 1. At the time the computations were made these 
values of the derivatives were considered those most representative of 
the configuration. For a condition with no spin there is no disturbance 
in the yaw plane and therefore the motion is confined to the pitch plane 
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only. The variations of dynamic pressure Q, of deceleration along the 
flight path n, of Mach number and of the envelope of the angle- 

of-attack oscillation a with altitude are shown in figure 4(a) . The 
variations of Q, n, and Nj^ are typical of those for most of the 

runs and therefore are not presented for the rest of the runs. The varia- 
tion of altitude with range is shown in figure 4(b); this variation also 
is typical of the trajectories for most of the cases. It should be men- 
tioned that since this is a plot of altitude against range the true 
curvature of the flight path in inertial space is not shown. In order 
to obtain the true curvature, a plot of x against y would be needed. 

At the beginning of the reentry trajectory, 380,000 feet, the dynamic 
pressure was so low that an aerodynamic moment sufficiently large to 
turn the vehicle along the curvature of the flight path was not generated, 
and consequently the angle of attack increased from the initial 1° up to 
1.6° at 374,000 feet. At this point the pitching moment became large 
enough to turn the vehicle along the flight path and an oscillation in 
angle of attack developed. The first few cycles of the oscillation are 
shown in the angle -of -attack envelope plot in figure 4(a) . The 
increasing dynamic pressure as the vehicle descended through the atmos- 
phere had a constraining effect on the amplitude of the oscillation. 

The fact that a minimum value of a occurs before the maximum value 
of Q occurs is due to the destabilizing effects of axial force and 
negative damping. Below the point of maximum Q the angle-of -attack 
amplitude expands rapidly until %a = !• Below %a = 1 the positive 

aerodynamic damping rapidly decreases the amplitude. 


Effect of Aerodynamic Derivatives on Amplitude Envelope 
of Angle of Attack Above a Mach Number of 1 


The effect of changes in the damping derivative the normal- 

force-curve slope stability parameter are 

shown in figure 5. Varying the damping coefficient from =0.37 

to ^m^ = - 0.90 (fig* 5(a)) decreased the amplitude of the angle-of- 


attack oscillation as would be expected. This effect is confined to 
the lower portion of the trajectory where the decreasing dynamic pres- 
sure permits the effect of the aerodynamic damping to develop. During 
the upper portion of the trajectory the confining effect of increasing 
dynamic pressure was predominant and little effect of aerodynamic 
damping is seen. A dynamic stability factor defined in reference 8 is 


given by K = 


-2C X + C Za + 



This parameter varied 


from 4.04 for C m ^ = 0.37 to -0.14 for C m ^ = - 0.90 (table I) 


00 \D t"— 
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indicating a change from strong instability to slight stability. It is 
interesting to note that only for the case where the dynamic stability 
factor is negative does the amplitude of the angle -of -attack oscillation 
continuously decrease from reentry down to the altitude where = 1. 


The damping contributed by the normal-force-curve slope C 7 is 

illustrated by the comparison shown in figure 5Cb)* For the case where 
the damping coefficient C m ^ was zero, it can be seen that a negative 

value of corresponding to a positive lift-curve slope has a sta- 

bilizing effect. 


The effect of changes in the static stability parameter Cn^ on the 


angle-of -attack envelope is shown in figure 5(c) and on the frequency of 
oscillation in figure 6. For the basic configuration the maximum fre- 
quency of oscillation was about 0.9 cycle per second at maximum dynamic 
pressure. The values of were altered sufficiently to decrease the 

maximum frequency to 1 /4 cycle per second in one case = “0.016^ and 

increase the frequency to 4 cycles per second in the other ^0^ = - 4.000 
For the lower frequency case (^m^ = -0.0l6) a considerably larger ampli- 


tude of oscillation was obtained throughout the trajectory; this was a 
result of the behavior at the early stage of the trajectory. Because 
of the considerably smaller value of tor this case the angle of 


attack built up to almost 3*5° compared with 1.6° for the basic condi- 
tion before the vehicle was able to develop a sufficient moment to follow 
the flight path. For the case in which the frequency of oscillation was 
higher the vehicle almost immediately starts to follow the flight path 
with the consequently smaller amplitude of oscillation throughout the 
run. The computations for this high-frequency case were stopped at 
about 170,000 feet because the high frequency of oscillation required 
that a very large number of data points be computed in order to define 
adequately the oscillation. This would consume an unusually great 
amount of computing time and the additional information to be obtained 
was not considered worth the machine effort required. In addition to 
the variation of frequency with altitude obtained from the complete com- 
putations, the frequency obtained from the simple relationship 


Frequency = 


i_ p c 4x 

2« V I 


is shown in figure 


6 . 


The agreement between the 


two methods is very good, which indicates that if the dynamic pressure 
is known the simple computations would give an adequate representation 
of the frequency of oscillation. 
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Effect of Hate of Spin 

The effect of constant rate of spin about the vehicle longitudinal 
axis on the angle -of -attack envelope is shown in figure 7. The motion 
of the vehicle for the spinning case is not an oscillation in the normal 
sense as it is for the nonspinning case but a precessional type of motion; 
the double curves in figure 7 indicate the maximum and minimum angles of 
attack during the precessing motion. The fact that the maximum and the 
minimum angles of attack are so nearly the same indicates that the pre- 
cession is nearly circular about the flight path. The data show that 
for a spin rate of 1 rpm, which is considered sufficient to n ull ify the 
dispersion due to configuration asymmetry, the maximum angle of attack, 
which occurs at the Mach number of 1, is only slightly larger (about 2°) 
than that obtained for the nonspinning case. Increasing the spin rate 
had a deleterious effect on the angle-of -attack amplitude in that for a 
spin rate of 19.1 rpm the angle-of -at tack amplitude exceeded 80° and for 
a spin rate of 57.5 rpm it exceeded 90°. These large angles of attack 
result mainly from the increased stiffness of the spinning body relative 
to the normal aerodynamic moment tending to turn the body into the flight 
path. Because of this increased resistance to turning and the low aero- 
dynamic moments during the early part of the trajectory, the spinning 
vehicle tended to remain in its initial attitude considerably longer than 
did the nonspinning vehicle; therefore, the spinning vehicle built up 
to a larger value of angle of attack. Thereafter, the effect of dynamic 
pressure and damping were similar to that obtained for the nonspinning 
case. The fact that the ratio of the maximum angle-of-attack amplitude 
obtained at Nf 4 a = 1 to the minimum angle-cf-attack amplitude is about 

the same for the spinning and nonspinning cases (ojnax/^tnin = 20 for 

the nonspinning case, 20 for the 1-rpm case, and 23 for the 19.1-rpm 
case) indicates that the instability obtained near Mach number 1 is 
attributable mainly to the factors discussed in the previous sections. 


Effect of Earth Spin 

The effects of including the geostrophic winds (winds due to the 
earth's spin) on the angle-of-attack envelope and trajectories are shown 
in figure 8. The data in figure 8(a) show tnat the geostrophic winds 
in the form of head or tail winds have a negligible effect on the sta- 
bility of the vehicle. For the crosswind case (vehicle traveling north, 
perpendicular to earth's spin) the angle-of-attack pattern is similar 
but approximately twice as large as it is for the other cases. In this 
instance the geostrophic wind adds an appreciable angle in the yaw plane 
to the angle already present in the pitch plane, whereas for the other 
cases the geostrophic wind merely alters the relative wind by about 

52 percent. There is, however, an appreciable effect on the location 
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of the impact point for all three cases. The trajectories in fig- 
ure 8(b) show that for a reentry in the same direction as the earth's 
s pi n (going east) there is a tail wind and the vehicle travels farther 
than for the nonspinning case (basic condition) . A point on the surface 
of the earth directly under the point of reentry will move in the same 
direction as the vehicle because of the earth's spin and thus the appar- 
ent range will be decreased. During the time the vehicle is descending, 
the point on the surface actually moves farther than the increase in 
range caused by the geographic tail wind, and, therefore, the net range 
is less for a body reentering in the same direction as the earth spin 
than it would be for a stationary earth. (See the following table.) 

For a body reentering against the earth's spin exactly the opposite is 
true . 


Direction of reentry 

Vehicle total 
travel, miles 

Travel of a 
point on surface 
of earth, miles 

Net vehicle 
travel, miles 

East (with earth spin) 

2,102 

206 

1,896 

No spin 

1,989 

0 

1,989 

West (against earth spin) 

1,884 
i _ 

202 

2,086 


For a vehicle heading north at reentry the basic trajectory is not 
appreciably affected by introducing earth spin (1,989 miles). The 
vehicle, however, drifts in an easterly direction 76 miles during its 
flight because of geostrophic winds. A point on the surface originally 
at the same longitude as the reentry point but at the latitude at which 
impact takes place will have moved eastward 180 miles. The net effect, 
therefore, is such that the vehicle appears to move westward 104 miles. 


Effect of Reentry Angle and Velocity 

The basic initial velocity for most of the computations was 
25,752 ft/sec which is the orbital velocity for a body in a 380,000- 
foot circular orbit. The data in figure 9 show the effect of increasing 
the angle of reentry from 0° to 3° for the same initial velocity and at 
a 1° angle of attack. The data in figure 10 show the effect of 
decreasing this initial velocity up to 10 percent for a reentry angle 
of 0°. The effect on the range traversed from initiation of the reentry 
to the point of impact was much the same for both variables in that 
there was a very rapid decrease in range for the initial 1° declination 
or the initial 2-percent reduction in velocity. Thereafter, the relative 
change in the range was much smaller. 




12 


Although all these runs started with an initial angle of attack of 
1° they exhibited an increase in angle of attack at the start of the 
trajectory because the dynamic pressure developed at this high altitude 
was insufficient to turn the vehicle along the flight path (figs. 9(c) 
and 10(b)). Once sufficient dynamic pressure has been developed, how- 
ever, the angle-of-attack envelope pattern is similar for all the runs, 
and the difference in magnitude depends upon the magnitude of the angle 
of attack attained at the beginning of the run. 


Effects of Single-Point-Attachment Drogue Parachute, of C ^ , 
and of on -Attack Envelope at Low Speeds 

The effect of the damping coefficient C mq for both the clean 

vehicle and the vehicle with a single-point-attachment drogue parachute 
is presented in figure 11. The beneficial effect of having good aero- 
dynamic damping is immediately apparent. If the damping is good 
^m^ = -0.5 or better ^ the drogue parachute helps very little. If, on 

the other hand, the damping is low the drogue parachute will contribute 
materially to reducing the amplitude of oscillation. Changing the value 
of c Z a fr °m a destabilizing value of 0.6 to a stabilizing value of -0.28 

produces an appreciable improvement when the damping is low and, as noted 
before, this type of change has only minor effect when good damping 
exists (fig. 12). 

In addition to this damping effect, the drogue parachute increases 
the frequency of oscillation by a factor of At an altitude of 
72,000 feet the frequency of oscillation increased from 0.20 cps with- 
out the drogue parachute to 0.6l cps with the drogue parachute. 


Effect of Two-Point Attachment of Irogue Parachute 

The effectiveness of the various methods of drogue -parachute attach- 
ment to the vehicle in decreasing the amplitude of the angle-of-attack 
oscillation for the case of C™ =0 is shown in figure 13 . The progres- 
sive benefits of going from no parachute to a single-point-attachment 
parachute then to the two-point-attachment of the parachute for 

^ = 0.3125 and for ^ = 0.625 are quite obvious. In the practical case 

the two-point-attachment methods mentioned would be at least four-point- 
attachment systems (two in the pitch plane and two in the yaw plane) . 

Each of the various parachute-attachment methods seems to act as a lim- 
iter on the amplitude of the oscillation. 
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Besides having an appreciable effect on the angle-of -attack ampli- 
tude the drogue parachutes have a marked effect on the frequency of the 
oscillation. This effect is shown in the following table where the fre- 
quencies for the various methods of parachute attachment are given for 
an altitude of 72,000 feet, which is just below the altitude at which 
they are initially employed; 


Type of attachment 

Frequency obtained 
from IBM data, cps 

Frequency obtained from 
one -degree-of -freedom 
computation, cps 

No parachute 

0.20 

0.21 

Single point 

.61 

.61 

Two point, a/d = 0.3125 

1.36 

1.36 

Two point, a/d = 0.625 

1.86 

1.85 


The increase in the frequency of oscillation as the configuration is 
changed from one with no parachute to one with an extended two-point 



attachment is obvious . 


The table also shows that the fre- 


quency can be adequately predicted by making a one-degree-of-freedom 
computation of the time history of the angle of attack using the Laplace 
transform method so that the initial conditions for the offcenter para- 
chute attachment can be included. 


The effect of using a two-point attachment ^ = 0 . 3125 ^ instead of 

a single-point attachment for several damping coefficients ^m^ = 

- 0 . 5 * and -l.o) is shown in figure 1 ^. It can be seen that for each 
damping coefficient the use of the two-point attachment instead of the 
single-point attachment resulted in an improvement in the angle-of - 
attack envelope. For high values of C m of course, the improvement 


is not nearly so great as it is for low values of damping. Computations 

for the two-point-attachment - = 0.625, cases with large values of 

d 

damping were not made because it was believed that the large improvement 
shown for C m =0 (fig- 13 ) indicated that there would be a progres- 
sive improvement for the larger damping values and that the machine time 
could be better utilized for other aspects of the problem. 
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Drogue -Parachute Effects at Large Angles of Attack 

In order to investigate the effects of the drogue parachute opening 
at angles of attack other than a value of 5° (already discussed) compu- 
tations were made for several cases in which the initial conditions were 
similar to those for the parachute opening at the original 5° angle of 
attack but with the angle of attack arbitrarily increased to 10°, 20°, 

^■0 ° > &nd 60°. It should be mentioned here that these computations are, 
strictly, not applicable for the large angles considered because the 
stability derivatives used were assumed to be constant with angle of 
attack and because, for the development of the drogue— parachute terms, 
the analysis was limited to small angles of attack. However, the results 
should be useful for indicating the gross effects of the drogue -parachute 
behavior. The data for these computations are presented in figure 15 j 
it can be seen that for zero damping coefficient the drogue parachute 
attached at a single point can only effect enough damping to maintain 
the initial angle of 60°. This result is in agreement with the data 
presented in figure 11 where it was shown that for zero damping coeffi- 
cient the drogue parachute permitted the angle of attack to build up to 
a limiting value of about 60°. For a higher damping coefficient of 
Cmq = -0.5 it can be seen that the angle-of-attack envelope was damped 

quite well and also that the damping was nonlinear, being greater at the 
large angles than it was at the smaller angles. In order to show this 
effect somewhat clearer, three curves are shewn for each of the initial 
angles of attack of 60° and 40° in figure 15(b). One curve is the curve 
obtained for the normal 5° amplitude case multiplied by factors of 8 or 
12 in order to give initial angles of 40° and 60°, respectively. The 
second curve is one in which the drogue -parachute term is omitted from 
that portion of the equations where it is multiplied by qd/2V R in the 

Y-moment equation and rd /2V R in the Z-moment equation. It was included 

in all the other parts of the equation. (See appendix A.) The third 
curve is the same as that presented in figure 15(a) which was obtained 
by using the full equations. It can be seen that the parachute is effec- 
tive in introducing damping even though those terms that would normally 
be associated with a damping coefficient have been omitted and that 
when these terms are included there is a considerable increase in damping. 
The difference between the expanded 5° curve and the normal computations 
is a measure of the nonlinearity of the damping at these higher angles. 

The effect of the drogue parachute attached at two points at these 
higher angles of attack is shown in figure 15(c) for an initial ampli- 
tude of 60° and in figure 15(d) for an initial amplitude of 40°. It can 
be seen that going progressively from a single-point attachment to a 

two-point attachment ^ = 0.J125 then | = 0.625^ made substantial 

increases in the damping of the angle-of-attack envelope. As a matter 
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of fact, for the — = 0.625 two-point attachment the damping was satis- 
d 

factory even for C m ^ = 0 (fig. 15 (c)). The data in figure 15 (d) for 

40 ° initial amplitude indicate a very similar pattern to that for 60° 
initial amplitude and will not be discussed further. 


CONCLUSIONS 


An analytical study of the dynamic behavior of a typical manned 
nonlifting reentry vehicle has indicated the following conclusions: 

1 . In the high-speed portion of the trajectory (above a Mach number 

of 1) a moderate amount of negative damping can be tolerated without too 
adverse an effect on the angle -of -attack oscillation. Below a Mach num- 
ber of 1, however, reasonably good damping is necessary - -0.5) in 

order to avoid divergence of the angle of attack. 

2 . The stability can be improved, at low speeds, by employing a 
drogue parachute. The effectiveness of the drogue parachute is pro- 
gressively increased as the attachment is changed from a single point 
on the longitudinal axis to several points around the periphery of the 
rear face and then to outrigger arms extending outward from the surface 
at the rear . 

3 * A moderate rate of spin, enough to remove dispersion caused by 
vehicle asymmetry (about 1 rpm), would not have a serious adverse effect 
on the vehicle stability. Higher rates of spin, however, would be 
deleterious . 

4 . Geostrophic winds did not have a serious effect on the vehicle 
stability; they did, however, appreciably alter the impact point. 

5 . Increasing the reentry angle or decreasing the initial velocity 
from basic orbital conditions had a similar effect on the range in that 
a small change in either the angle (l°) or velocity (2 percent) at first 
effected a large decrease in range but thereafter similar changes had 
considerably smaller effects. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field., Va., April 11 , i960. 
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APPENDIX A 
EQUATIONS OF MOTION 


The equations of motion used in this investigation are for a modi- 
fied body axes system (fig. 1). The axes system differs from the usual 
body axes system in that the Z-axis is constrained to remain in a plane 
parallel to the K^Z^ inertial plane and in that the body is free to 
spin about the X body axis. The X-axis is always alined with the axis 
of symmetry of the body. These equations are in a form frequently 
employed in ballistics work, and they may be derived by resolving equa- 
tions of motion in a normal body axis system into the XYZ modified body 
axis system as shown in reference 11. Therefore, a full development of 
the equations will not be presented. The equations used herein are as 
follows : 


X-force equation: 


u ■ #( c x' * c XjP ’) 


+ rv - qw 


R e . 

-g — T— (x cos \|r cos 0 + y sin 0 + z sin i|r cos 0) 
R* 


where Cy' = Cy cos a and ^x,p' = Cy p coe a - 
Y-force equation: 


* ' wl( c r a + °x,p)( T ’ V 


+ \Cv - 4c 


*R 


x,p 


' . 2 

l/ v ~ V Y) , a* A 

a \ v R ) d { 


V, 


rd 

2Vr/ 


R 

- ru + qw tan 0 - g -|-(-x cos \|r sin 0 + y cos 0 - z sin \|r sin 0) 

R^ 


Z-force equation: 


* ■ ir(( c z a + c x,p >)(^T/ + + kc *,i 


R, 


l( w ~ V z\ 2 . a**/ w “ v z\ 
d \ V R / d \ V R j 


qd_ 

2V R J 


+ qu - qv tan 0 - g (-x sin + z cos \|f) 

R 5 


1-q 00 VO t>- 



•^1 00 VO o- 


3T 


Y-moment equation: 



I X 

- — pr + rq tan 0 
I Y 


Z-moment equation: 



^ = p - q tan 0 
0 = r 


^ = -q/cos 0 
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and 

x = u cos 0 cos + - v sin 0 cos - w sin \|r 

y = u sin 0 + v cos 0 

z = u cos 0 sin y - v sin 0 sin ^ + w cos \|r 


V Y - HRp cos 0 cos 


- ifjcos i cos 0 - sin^- “ tjsin ^ 


Vy = ftR e cos 0 


V z = HR e cos 8 


0 = arc tan 


- cos^ “ ^jcos \| f sin 0 + sin^ - \jrj 
(2 "" ^ sin ^ “ sln ^ " ♦jcoB \|rj 


sin 


- cos 


& 


+ z c 


= arc tan — 
z 


(!'*)' 

Range in equatorial plane = R g ^ - 

Range in meridian plane = R (§) 


cos 9 


i|r sin 



^ CO VO 
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APPENDIX B 


DEVELOPMENT OF TERMS USED TO SIMULATE DROGUE 
PARACHUTE IN EQUATIONS OF MOTION 


The terms used to simulate the drogue parachute in the equations 
of motion are developed in this appendix. In order to simplify the 
equations it is assumed that the line attaching the parachute to the 
vehicle is long relative to the magnitude of the displacement of the 
attachment point; therefore the drag force acts straight back in line 
with the relative wind, and any motion of the parachute is omitted. In 
addition, higher order terms which would considerably complicate the 
computations but add little to the effects investigated are omitted. 

The drogue parachute thus provides pure drag force and the forces 
and moments along the body axes are provided by the angular displace- 
ment of the body and the motion of the attachment point about the cen- 
ter of gravity of the body. 

The force applied to the vehicle by the drogue parachute is given 

*>y 


F P " VS? C DjP 

and the moment is given by 

M P = % A p C D,p k p 

where kp is the moment arm of the parachute force. The dynamic pres- 
sure Q_ acting on the parachute is 


^ - K 2 


By referring to the following sketch 
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Z 



it can be seen that 


- + AV = Vp + Z0 sin 0 f a0 cos 0 


and 


v p 2 = v r 2 + 2V R (z3 sin 0 + a§ cos 0) + Z 2 0 2 sin 2 0 

• 2 .2 

+ 2al0 sin 0 cos 0 + a 2 0 cos 2 0 

Assume angles small enough so that sin 0^0 and cos 0 « Z. The equa 
tion for V p 2 thus reduces to 

v p 2 = v r 2 + 2V r 0(Z0 + a) + 0 2 (z 2 0- + 2aZ0 + a 2 ) 


The dynamic pressure becomes 


S 




i + 4 



+ 4 


• 2 2 / 2 2 
£JL( 6 + 

4V R 2 d 2 



The term 



is a small fraction and therefore it is believed that 


omission of the terms containing 



would not affect the results 


hi 00 vo t- 



-3 G\ cd \r* 
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appreciably. Thus, the equation for the dynamic pressure can be written 
as 

% ■ H 2 [a ♦ * e ♦ |j 

The tern kp of the moment equation is k p = l sin 0 + a cos 9 = 29 + a. 
The equation for the force now becomes 

■ • 

F_ = QA-^C-p. 1 + 4. ®!-[i 9 + 

P D,P[ 2V R Vd d ) 

and the equation for the moment becomes 


Mp - QApCj^p Id + a + 4 


The forces and moments due to the off center attachment points of 
the drogue -parachute lines (the terms containing a) require some further 
explanation. As the vehicle oscillates from positive to negative angles 
of attack first one parachute line will be taut and the other slack; 
then, they will reverse. The taut line bears the full load whereas the 
slack line has no load. Therefore, at the angle of attack where the 
lines change from slack to taut there will be an immediate shift in load 
from one side to the other and the load will not be proportional to the 
angle of attack. Thus, it can be recognized that this term is a step 
function which changes sign with the angle of attack. 

The force and moment coefficients to be used in the equations can 
now be determined. 

X-force coefficients: For this force the second and third terms 

in the force equations ^those containing the term are small com- 

pared with 1 and can therefore be neglected and the equation for the 
incremental force along the X-axis due to the parachute reduces to 



AF X = QA AC X = -QApC D>p 
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and for ACy the following expression is ottained: 


ACv = - 




QA 


Let — C n ^ 
A D >P 


= -C 


x,p- 


Thus , 


= c X,p 


Y-force coefficients: The equation for the incremental force along 

the Y-axis due to the parachute is 


AFy - QA AC y - QApC D#p jl + 4 6 + || 


and 


ACy — 


-c x [e + 4 q 2 + S: ej] 

X,P L 2Vp\d d /j 


The coefficient can be broken up into two terms, a term ACy propor- 

x <x 

tional to a and a term ACy proportional to . Taking into 

r 2 V r 

account the proper signs for the axes system used gives the following 
expression for ACy 
x a 


AC Y a = c X,p 


and the term ACy^ becomes 


AC Y = -4Cy^ p 


i fLi h 

dV V R 


ft 


-v) G\ CD fcH 
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Z-force coefficients: In a manner similar to that used to obtain 

ACy it can be shown that the terms in the Z-force equation are 


= C X,p 



Y-moment coefficients: The incremental moment about the Y-axis 

due to the parachute is 


£M y = -QAd AC m = -QApC^p 


~i 0 + 8 . + 


9d (l‘ 
2V 0 \d 


j2' 

9^ + 2 — 9 + — , 
2V^Vd d d 


,2 . ^ aZ „ . a c 


and the incremental coefficient becomes 


^m = C X,p 


i e + £ + 4 e 2 + 2 2i e + £ 

a d ^Ma 2 a 2 a 2 / 


Here, as for the forces, the total coefficient can be broken up into 
sections so that when proper cognizance of signs is made the coeffi- 
cients become 


AC 


% 



Z_ 

d 


and 




2k 


Z-moment coefficient.: In a manner similar to that used to obtain 

the Y-moment coefficients it can be shown that 


AC 


% ‘ 



d 



AC n 


_r a 

%p d 


CO VO 



1^ CO VO t- 
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TABLE II . - MASS AND DIMENSIONAL DATA 


m. 


“X* 


Y' 

A, 

l/d 

a/d 


slugs • 
slug-ft 

slug-f t 

sq ft . 

ft • . 


2 

2 


W 


(above M = 1 ) 


Condition I 

Condition II 

53-2 

65.22 

448 

280 

473 

440 

34.91 

34.91 

6.667 

6.667 

1.0125 

1.0125 
0.5125 and 0.6250 

31.66 

38.81 
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Figure 4.- Concluded. 
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Figure 5*- Effect of aerodynamic derivatives on angle-of -attack envelope. 
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Figure 5.- Continued. 
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Figure Concluded. 
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Figure 6.- Variation of computed frequency of oscillation with altitude. 



Rate of spin, rpm 


39 



minimum angles of attack. 


Basic condition no earth spin 
Reentry in same direction as earth spin 
Reentry in opposite direction to earth spin 
Reentry perpendicular to the earth spin 
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Figure 8.- Effect of geostrophic wind on angle-of -attack envelope and trajectory. 
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(b) Effect of geostrophic wind on range for head or tail winds. 

Figure 8.- Continued. 
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Figure 8.- Concluded. 
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(a) Variation of range with angle of declination of reentry. 


Figure 9.- Effect of initial declination of flight path on range trajec- 
tory and angle-of -attack envelope for orbital initial velocities. 
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Figure 9*- Continued. 
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(a) Variation of range with initial velocity. 


Figure 10.- Effect of change in initial velocity on range and angle-of- 
attack envelope for an entry angle of 0°. a^, = 1°. 
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(b) Effect of initial velocity on angle-of-attack envelope. 
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Figure 11.- Effect of on angle -of -attack envelope below N^ a = 1 without and with single- 

point-attachment drogue parachute. 
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Figure 13- Comparison of angle -of -attack envelopes for two-point and single -point parachute 
attachments with the envelope obtained when no drogue parachute is attached. C m =0. 
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Figure 15 •- Continued 
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Figure 15 •- Continued. 
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Figure 15 ■- Concluded. 



